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Abstract 

Jet production in charged and neutral current events in the kinematic range of Q 2 from 
640 to 35 000 GeV 2 is studied in deep-inelastic positron-proton scattering at HERA. The 
measured rate of multi-jet events and distributions of jet polar angle, transverse energy, 
dijet mass, and other dijet variables are presented. Using parton densities derived from 
inclusive DIS cross sections, perturbative QCD calculations in NLO are found to give a 
consistent description of both the neutral and charged current dijet production. A direct, 
model independent comparison of the jet distributions in charged and neutral current events 
confirms that the QCD dynamics of the hadronic final state is independent of the underlying 
electroweak scattering process. 
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1 Introduction 



Deep-inelastic scattering (DIS) at the electron-proton collider HERA offers unique possibilities 
to reveal the partonic structure of matter. At very high four-momentum transfer squared — Q 2 
the exchange of all the electroweak gauge bosons (photon, Z° and W /± ) becomes important 
allowing the standard model of electroweak and strong interactions to be tested at distances 
as small as 10 _18 m. The inclusive DIS cross sections of neutral current (NC) ep — > eX and 
charged current (CC) ep — > vX interactions have been measured [[]], and are well described 
by the standard model. In this analysis we complement these results by the first detailed inves- 
tigation of dijet structures in both NC and CC processes. 

Within the Quark-Parton-Model DIS gives rise to events with (1+1) jets, where one jet 
originates from a quark struck out of the proton and a second jet is due to the proton remnant 
(denoted '+1'). Events with (2+1) jets, referred to as dijet events, are predicted by Quan- 
tum Chromodynamics (QCD) due to contributions in 0(a s ), namely QCD-Compton scattering 
eq — > eqg and Boson-Gluon-Fusion eg — > eqq as illustrated in Figure []]. In CC interactions 
several events with multijet structures have been identified [J3]] and the jet shape has been mea- 
sured [Q]. However, due to the relatively small number of CC events observed so far at HERA, 
the structure of the hadronic final state has not yet been studied in detail. In NC interactions 
clear multi-jet structures have been established [|5|] and have been used to test QCD [§]. Pre- 
vious analyses of dijet production in NC processes, however, did not yet extend to very high 
values of Q 2 . 

In the present paper a dijet analysis of a sample of 460 CC events and approximately 8 600 
NC events with Q 2 in the range of 640 to 35 000 GeV 2 is performed. Various dijet distributions 
are compared with the predictions of QCD Monte Carlo models and with perturbative QCD 
calculations in next-to-leading order (NLO). In addition, the jet distributions of the CC and the 
NC events are compared directly, in order to test the hypothesis that QCD radiation proceeds 
independently of the underlying electroweak scattering process. 

2 Detector description and data selection 

This analysis is based on the data sample recorded with the HI detector in the data taking 
periods 1994-1997 at HERA. In this period HERA was operated with positron and proton 
beams of 27.5 and 820 GeV respectively, corresponding to a centre-of-mass energy of ^fs « 
300 GeV. The collected integrated luminosity for this analysis is 35.6 pb _1 . 

2.1 Detector and trigger 

The components of the HI detector 0] most relevant for this analysis are the central tracking 
system, the liquid argon calorimeter and the instrumented iron return yoke. 

The central tracking system consists of two concentric drift chambers covering a polar an- 
gular rangeQ of 15° to 165°. Two polygonal drift chambers with wires perpendicular to the beam 

'The forward direction and the positive z-axis are defined as the proton beam direction. The origin of coordi- 
nates is the nominal ep interaction point. 
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axis improve the determination of the z coordinate of the measured tracks. The central tracking 
system is surrounded by a liquid argon sampling calorimeter covering a polar angle range of 
4° < 9 < 154°. The electromagnetic and hadronic sections of the liquid argon calorimeter 
correspond in total to a depth of 4.5 to 8 interaction lengths. The energy resolution of the liq- 
uid argon calorimeter for electrons and hadrons was determined in test beam measurements to 
be a/E = 12%/ v /E(GeV) © 1% and cr/E = 50%/ v /E(GeV) © 2%, respectively [§]. The 
systematic uncertainty of the electromagnetic energy scale is determined to be 0.7% for the 
majority of the selected events and increases to 3% at the highest Q 2 [p. The uncertainty on 
the hadronic energy scale of the liquid argon calorimeter is 4%. 

Outside the calorimeters a large superconducting solenoid provides a magnetic field of 
1.15 Tesla. The instrumented iron return yoke identifies energetic muons and detects leakage of 
hadronic showers. 

The trigger conditions for CC events are based on the reconstruction of a large missing 
transverse momentum in the trigger sums of the liquid argon calorimeter NC events are 
triggered on the basis of a localized high energy deposit in the electromagnetic part of the 
calorimeter [|7|]. 

2.2 Event selection 

Selection of CC events 

The selection of CC events is similar to those of [[I], |TOp. It is based on the observation of a large 
imbalance in transverse momentum due to the antineutrino escaping direct detection. The trans- 
verse momentum P^ ad , reconstructed with the liquid argon calorimeter and the instrumented 
iron, is required to exceed 25 GeV. No scattered positron must be found in order to reject neu- 
tral current events. The z coordinate of the primary event vertex z vtx has to be within a distance 
of ±35 cm from the nominal ep collision point. The inelasticity yhad = Eh(l — cos 9h)/2E e , 
calculated from the energy depositions in the calorimeters and the energy of the positron beam 
E e , must be in the range 0.03 < yhad < 0.85. The kinematic selection criteria imply a minimum 
virtuality Q 2 of the exchanged boson of 640 GeV 2 . 

Background events due to cosmic muons, beam-halo muons and beam-gas interactions are 
removed by further requirements on the event topology and timing JlOt]. Furthermore a visual 
scan of the remaining events is performed. 

The final event sample consists of 460 CC events. The background from photoproduction 
events is less than 2%. It is estimated from Monte Carlo simulations and from data events where 
the scattered electron is detected at very small scattering angles. The number of background 
events from other sources is negligible. 

The trigger efficiency for events with P£ ad > 25 GeV has been determined as a function of 
the kinematic variables and the jet variables studied using NC events where the information of 
the scattered positron [ |10| ] is discarded. The average trigger efficiency is « 95%. It is corrected 
for in all measured distributions. 

Selection of NC events 

The kinematic selection criteria of the NC events correspond to those of the CC events. The 
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NC selection requires the identification of the scattered positron. Fiducial cuts are applied to 
the impact position of the scattered positron in the liquid argon calorimeter in order to avoid 
inhomogeneities at the boundary of detector modules [|]|. The kinematic selection is based on 
the variables pj, and y e reconstructed from the scattered positron momentum with the exception 
of the requirement yhad > 0.03. The summed energy E and longitudinal momentum compo- 
nents P z of all reconstructed detector objects (see section 4.1) must fulfil E — P z > 35 GeV 
to suppress QED radiative events. The z coordinate of the primary event vertex is required to 
be in the same range as for the CC events. The NC sample consists of ~ 8 600 events with a 
negligible number of background events from photoproduction. The trigger efficiency of the 
NC events is ^99%. 

The distributions of the reconstructed kinematic variables P T , Q 2 , the Bjorken scaling vari- 
able x and y for the CC and the NC event samples are shown in Figure 0. A good description of 
the data by the Monte Carlo model ARIADNE (see next section) combined with the HI detector 
simulation is observed. Note that the distributions are normalized to the total number of CC or 
NC events Nnis respectively. The differences between the CC and NC distributions [jT]] are due 
to the different couplings and propagators of the bosons in CC and NC interactions. 

The main selection criteria for the CC and NC events are summarized in Table [[]. 



CC 


NC 


no e + found 

P£ ad > 25 GeV 
0.03 < y had < 0.85 


e + found 

V e T > 25 GeV 
0.03 < y had , y e < 0.85 


\z vtx \ < 35 cm 


E - P z > 35 GeV 

\z vtx \ < 35 cm 



Table 1 : Selection criteria for the CC and NC DIS event samples 



3 QCD Monte Carlo models and QCD NLO programs 
3.1 QCD Monte Carlo models 

Four different QCD Monte Carlo models are used in this analysis: ARIADNE 4.10 Q; HER- 
WIG 5.9 Hg], LEPTO 6.5.2/5 [fj and RAPGAP 2.08/06 0. All models use the LO matrix 
elements for QCD-Compton and Boson-Gluon-Fusion. ARIADNE implements higher order 
QCD processes with radiating colour dipoles |JT3|], HERWIG, LEPTO and RAPGAP use initial 
and final state parton showers instead [|T6|]. In the context of this analysis LEPTO and RAPGAP 
are similar. They differ in the way the divergences of the LO matrix element are regulated. Frag- 
mentation of partons into hadrons is modelled with the Lund string model [[FT]] in ARIADNE, 



LEPTO and RAPGAP, and with the cluster model [18] in HERWIG. 



The latest versions of the models as described in [[TPI] are used. The LEPTO version used 



contains a refinement of soft-colour interactions, the generalised area law model [20]. LEPTO 



6 



has been tuned to describe jet distributions at HERA and the corresponding values of the model 
parameters are taken here. In HERWIG, we use the leading order and not the next-to-leading 
order formula for a s as proposed in [JT5]. The parton density functions CTEQ4L pT| ] are taken. 

ARIADNE and LEPTO are incorporated into DJANGO [[22p, version 6.2, which simulates 
the effects of QED radiation. 



3.2 QCD NLO programs 

Four programs MEPJET |JZ3p, DISENT [0J], DISASTER++ @] and JETVIP are avail- 
able for perturbative QCD calculations of jet cross sections in NLO. Currently MEPJET is the 
only NLO program that considers W or Z° exchange. MEPJET is thus used to calculate the jet 
distributions in CC processes. The NC jet cross sections are calculated with DISENT follow- 
ing the recommendations in fl27p. We use the parton density functions determined by the HI 
Collaboration [[p and choose Q 2 as the renormalization and factorization scales, n 2 R and (j, F . 

We compared the predictions of MEPJET and DISENT for the jet distributions presented 
below. In leading order, we find agreement within a fraction of a per cent. In NLO, MEPJET 



is systematically lower than DISENT by ps 10% confirming the results of [[281] . Note that the 
comparison of various jet cross sections calculated with DISENT and DISASTER++ showed 
good agreement [EBfl. The differences between DISENT and DISASTER++ observed for ex- 



treme values of event shape variables []29|, |30|] are not relevant to this analysis. Currently it is 
unknown if the inconsistencies between the NLO programs observed in NC influence the CC 
predictions as well. 

DISENT does not consider Z° exchange, which for Q 2 > 5 000 GeV 2 reduces the inclu- 
sive e + p NC cross section by less than 5% on average compared with purely electromagnetic 
exchange. Since the dijet cross sections are also reduced, the effect on jet distributions normal- 
ized to the number of DIS events is small. Correction factors were calculated using ARIADNE 
4.10 and are applied to any DISENT prediction at Q 2 > 5 000 GeV 2 . 

In order to compare the perturbative QCD predictions to the data, bin-by-bin hadroniza- 
tion corrections are determined using the QCD models ARIADNE and HERWIG. The average 
correction factors from the two models are applied to the NLO distributions. The maximum 
deviation between the average correction factor and the correction factor for either model alone 
is taken as hadronization uncertainty. 



4 Definition of jet observables 

4.1 Jet algorithm 

Jets are reconstructed with a modified version of the Durham jet algorithm which was originally 



introduced in e + e annihilation experiments [3T]. The algorithm is applied in the laboratory 



frame. It is modified for application in DIS in two respects: a missing-momentum four-vector 
is introduced which is treated as an additional object by the jet algorithm to account for the 
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momentum carried by the proton remnant escaping through the beam pipe; in NC events the 
scattered positron is removed from the final state objects and is only used to determine the 
missing-momentum vector. In CC events this is achieved by first reconstructing the neutrino 
from the hadronic final state, exploiting energy and momentum conservation. 

The jet algorithm calculates the quantity k\ ^ = 2iain[Ef, E 2 ] (1 — cosfly) of pairs of 
objects or 'proto' jets i, j. Here E^ and Ej are the energies of the objects i and j, and is 
the angle between them. The jet algorithm combines the pair of objects i, j with the minimum 
k\ f . to be a 'proto' jet by adding their four-momenta pi and pj. This prescription is repeated 
iteratively for the remaining objects until exactly (2+1) jets remain. At this stage, every event is 
treated as a dijet event by definition. Pronounced dijet structures are then selected by imposing a 
lower limit on y 2 , defined as the minimum k\^jW 2 of any combination of the (2+1) jets. Here, 
W is the invariant mass of the hadronic final state. It is calculated from all objects entering the 
jet algorithm. 

In order to determine the fraction of events with say (1+1) or (3+1) like jet structures it is 
more convenient to run the algorithm with a fixed jet resolution parameter y cut . In this case the 
iterative clustering procedure ends, when the k\^jW 2 of any pair of objects or proto jets is 
larger than a given value y cut . Thus the number of jets reconstructed varies from event to event. 

The algorithm is applied to the tracks reconstructed in the central tracking chambers and 
the energy depositions (clusters) in the liquid argon calorimeter. For tracks and calorimeter 
clusters that can be matched, the energy is determined from either the calorimeter alone or from 
a combination of track and cluster energy as described in The polar angle of each detector 
object, track or cluster, is required to exceed 7° in order to avoid the region close to the edge of 
the calorimeter. This improves the resolution of the reconstructed jet quantities. 

The same definitions of the jets are used for the analysis of the data and the Monte Carlo 
events after detector simulation. In events simulated at the hadron or parton level and in the 
perturbative QCD calculations, the jet algorithm is applied to hadron or parton four-momenta 
respectively. The polar angle cut of 7°, which is applied for detector objects, is also applied for 
hadrons and partons. 

In the calculation of y 2 the effects of the hadronic energy scale uncertainty largely cancel due 
to the method chosen to reconstruct W. The choice of a jet algorithm working in the laboratory 
frame leads to reduced experimental errors since a boost into another frame is avoided. This is 
relevant for the CC events where the resolution of the kinematic variables is worse than in NC 
events. 

In the present analysis the jet polar angles must fulfil 10° < 9j et < 140°. This restricts the 
jets to the acceptance of the liquid argon calorimeter. 

4.2 Jet observables 

The rate of events with jet multiplicity i is defined as Ri(Q 2 ) = Ni(Q 2 )/N DIS (Q 2 ), where 
Ni=i ; 2,3 is the number of events with one, two or three jets, not counting the proton remnant jet. 
Ndis is the number of selected deep-inelastic events. The value of y cut is taken as 0.002. 
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The dijet sample is correspondingly defined by the requirement y 2 > 0.002. Large values 
of y 2 correspond to events with (2+1) jets that are clearly separated and indicative of hard QCD 
radiation. Small values of y 2 are typical for events which intuitively may be considered as 
(1+1) jet events. In Figure [3] two CC events with very different values of y 2 are displayed for 
illustration. Note that, with this definition of the dijet sample, the events contributing to the 
one-jet rate as introduced above are eliminated, and the few three-jet events are now treated as 
dijet events. 

The distributions of the dijet variables y 2 , m 12 , z p , x p , the polar angle 9f wd and the transverse 
energy Exjwd of the most forward (non-remnant) jet are studied. m\ 2 is the invariant mass of 
the two non-remnant jets. The variables z p and x p are defined by 



where and 9i are the energies and polar angles of the two (non-remnant) jets remaining after 
the clustering of the jet algorithm. The variable z p corresponds to 1/2 min(l — cos 9*) where 

( 2=1,2 

9* is the polar angle of the parton i in the centre-of-mass system of the virtual boson and the 
incoming parton. In leading order QCD x p is equal to the ratio x/£ where £ is the fraction 
of the proton's four momentum carried by the parton entering in the hard scattering process 
(see Figure [T]). In the limit where one jet is absorbed into the remnant jet z p approaches 0. In 
the other limit, where the two (non-remnant) jets become one jet, m\ 2 approaches and x p 
approaches 1. 



5 Correction of the data 

The data are corrected for the effects of detector acceptance and resolution, and of QED ra- 
diation using the QCD models ARIADNE and LEPTO. For the correction of CC events, the 
number of events simulated for either model is approximately 150 times larger than that of the 
experimental data. The number of simulated NC events is at least six times larger than that of 
the data. The same event cuts and track/cluster selection criteria are applied to the simulated 
events and to the data. 

Correction of detector effects 

The measured jet distributions are corrected for detector effects with bin-by-bin correction fac- 
tors. The purity, defined as the number of simulated events which originate in a bin and are 
reconstructed in it, normalized by the number of reconstructed events in that bin, is on average 
60% for both CC and NC distributions. The purities estimated with LEPTO and ARIADNE are 
very similar. 

The stability of the results was tested by correcting the jet distributions reconstructed from 
simulated LEPTO events with the correction factors derived with ARIADNE. The agreement 
of the corrected jet distributions with the 'true' LEPTO jet distributions is good. Deviations are 
typically of a few percent. The largest deviations of 10 — 15% are seen in the z p distribution. 
These effects are considered in the model uncertainty discussed below. 




Q 2 



Q 2 + m 2 v 
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Correction of QED radiative effects 

The effects of QED radiation are considered by correction factors also. These factors are ob- 
tained from the ratio of the Monte Carlo distributions generated with and without inclusion of 
QED effects. The size of the corrections is ~ 5% for both CC and NC distributions. 

The combined detector and QED radiation correction factors from LEPTO and ARIADNE 
are averaged, and the resulting mean correction factors are used to correct the jet distributions. 

6 Determination of systematic errors 

The major sources of systematic errors are the model dependence of the detector corrections 
and the uncertainties of the electromagnetic and the hadronic energy scales of the liquid argon 
calorimeter. The total error of the majority of the CC data points is dominated by the statistical 
errors whereas the statistical and systematic errors are roughly of the same size for the NC data. 

Model dependence of correction factors 

The difference between the (average) corrected distributions and the distributions corrected with 
either model alone is taken as the error. The error is on average ~ ±3% for both the CC and 
NC distributions. 

Electron energy calibration 

The energy scale of electrons measured in the liquid argon calorimeter is known to 0.7% in the 
angular region where most events are situated. The effect on the jet distributions of the NC 
event sample is generally smaller than one per cent. 

Hadronic energy calibration 

In order to estimate the effect of the hadronic energy scale uncertainty on the measured jet 
distributions, the analysis is repeated with the hadronic cluster energies shifted by ±4%. The 
size of the corresponding changes depends considerably on the observable studied. For the 
CC events, the largest variation of ~ 11% is observed for the mvz distribution, the smallest 
variations of ~ 3% are observed for the y 2 , z v and 9f w d distributions. Similar variations are 
observed in the NC jet distributions. A variation of the track momentum by ±3% has negligible 
effect on the jet distributions. 

7 Results 

7.1 Jet event rates 

The rates of events with one, two and three jets Ri(Q 2 ) are shown in Figure as a function 
of Q 2 and are listed in Table |[ The jet event rates for CC and NC events are similar. For the 
chosen jet resolution, the fraction of one-jet events is « 70% for both CC or NC events. The 
fraction of dijet events is « 20% and that of three-jet events is a few percent. No CC three-jet 
events are yet observed at Q 2 > 5000 GeV 2 , which is statistically consistent with the QCD 
model expectations. The Q 2 dependence of the jet event rates is small. Note that Ri has a weak 
Q 2 dependence since most of the DIS events are reconstructed as (1+1) jet events. The jet event 
rates are well described by the QCD model ARIADNE. 



10 



7.2 Differential dijet distributions 



The measured CC dijet distributions of y 2 , mi 2 , z p , x p , Et, fwd an d 6f w d, corrected for detector 
effects and the effects of QED radiation, are shown in Figure |[ The distributions are based 
on the 120 CC events that pass the requirements y 2 > 0.02 and the 10° < 6 fwd < 140°. The 
measured differential dijet cross sections e.g. da dijet / d]j2 are normalized by the inclusive DIS 
cross section cr DIS for the kinematic selection of section 2. The y 2 and m 12 distributions are 
steeply falling. In the tails of these distributions events with clear dijet structures (see Figure |3|) 
and with dijet masses up to ~ 100 GeV are observed. The z p distribution shows a drop in the 
first bin at small z p which is due to the jet selection cut. The x p distribution is strongly peaked 
at large values of x p because the minimum Q 2 of the event selection is large. The corresponding 
average value of £ is ~ 0. 1. The forward jet distributions are strongly increasing at small polar 
angles and small transverse energies as is qualitatively expected by gluon bremsstrahlung off an 
incoming quark. 

In Figure ^] the corresponding distributions are shown for NC events. Here « 1900 events 
remain after the jet selection. The NC distributions show the same features as the CC distribu- 
tions. Note that due to the reduced statistical error of the NC distributions their total error is 
much smaller than that of the CC distributions. The differential CC and NC dijet distributions 
presented here are listed in Tables |3|. 

The corrected jet distributions are compared with the QCD models ARIADNE, HERWIG, 
LEPTO and RAPGAR Within the errors, the data are reasonably well described by the QCD 
models ARIADNE, HERWIG and RAPGAP. LEPTO roughly follows the data distributions but 
overall it is inferior to the other models. Significant deviations from the data are observed in the 
z p distribution in particular. These observations are valid for both CC and NC distributions. 

7.3 Comparison with perturbative QCD calculations in NLO 

The differential dijet distributions are also compared with QCD predictions in NLO. Two 
sources of theoretical error on the QCD predictions have been considered: the uncertainty of 
the hadronization corrections and the renormalization scale uncertainties. The size of these un- 
certainties is similar for CC and NC events. The hadronization corrections are typically smaller 
than 10%. Their uncertainty is estimated by the spread of the predictions of ARIADNE and 
HERWIG. The renormalization scale uncertainty of the NLO prediction is estimated by varying 
the renormalization scale ji 2 R from Q 2 to 1/4 Q 2 and 4 Q 2 . The resulting uncertainty is 5%. 
The non-remnant jets' average transverse energy in the Breit frame (E^ reit ) is ~ 10 GeV. 
Choosing (E^ reit ) 2 as renormalization scale changes the NLO predictions by « 5% compared 
with the scale Q 2 . A variation of the factorization scale has a marginal effect. 

The CC jet distributions are compared with the NLO calculations of MEPJET in Figure |7|. 
The MEPJET predictions provide a reasonable description of the data within errors. The cor- 
responding NC distributions shown in Figure |8] are well described by the NLO predictions of 
DISENT. 
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In the Figures [7] and [8| the NLO predictions for quark- and gluon-induced processes are 
also shown separately!. The predicted fraction of gluon-induced dijet events is ~ 20% in CC 
and « 15% in NC scattering for the selection criteria of this analysis. These fractions change 
by less than one per cent when varying the factorization or renormalization scale in the QCD 
calculations from Q 2 to 1/4 and 4 Q 2 . The dominance of quark-induced processes is mostly due 
to the relatively large values of x covered in this analysis. Figure [7] suggests that both quark and 
gluon contributions - calculated using the parton density functions determined from inclusive 
measurements - are needed to give a consistent description of dijet production in CC processes. 

The CC and NC jet distributions have also been determined for the subsample of the selected 
DIS events with Q 2 > 5 000 GeV 2 . After the jet cuts 17 CC and 91 NC dijet events remain. 
In Figure ^] the corresponding dijet mass distributions are compared with the NLO predictions 
of MEPJET and DISENT. Agreement is found in all bins. The measured distributions are also 
listed in Table |. 



7.4 Direct comparison of CC and NC dijet distributions 

In Figure [TO] the jet distributions of the CC events (full circles) are compared with those of 



the NC events selected in the same kinematic range (histogram). Systematic differences be- 
tween the jet distributions are observed in several bins. This is expected due to the different 
electroweak couplings and gauge boson propagators which also lead to different kinematic dis- 
tributions (see Figure Q). In order to account for these effects and to make possible a direct, 
model independent comparison of jet production in CC and NC processes, the NC events are 
reweighted to match the x and Q 2 spectra of the CC events. The weights are given by the ratio 
of the inclusive NC and CC DIS cross sections at the x and Q 2 of the NC event considered. The 
cross sections are calculated with DJANGO, and QED radiative corrections are taken into ac- 
count. Note that the inclusive cross sections do not depend on the hadronic final state properties 
and thus the reweighting procedure is independent of the modelling of the hadronic final state. 

The effect of this procedure has been tested with the QCD models ARIADNE, HERWIG 
and RAPGAR and the jet distributions of CC and reweighted NC events are predicted to agree 
within a few per cent typically. Residual differences between CC and reweighted NC jet dis- 
tributions are expected due to the different fraction of gluon induced events in CC and NC 
processes, to helicity effects [03L [34"|], and to the different parton densities contributing to NC 



and CC scattering. The NC event selection was repeated with the cut P^ ad > 25 GeV instead 
of p e T > 25 GeV. The changes in the corrected jet distributions are of the order of 2%. 

The measured NC jet distributions after reweighting are also shown in Figure [IT]. They are 
found to be consistent with the CC distributions confirming that at short distances the QCD 
dynamics of the hadronic final state are essentially independent of the underlying electroweak 
scattering process as is expected within the standard model. 

2 Note that the normalizations (Jdis of the quark-induced, gluon-induced and combined dijet cross sections are 
identical. 
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8 Summary 



A sample of 460 CC and of ps 8 600 NC events produced in deep-inelastic e + p scattering at 
HERA with the boson virtuality ranging from approximately 640 < Q 2 < 35 000 GeV 2 has 
been selected. In this sample jets are reconstructed using a modified version of the Durham 
algorithm. Jet studies are hence extended into a kinematic region where charged and neutral 
gauge bosons contribute at comparable level. 

Events with dijet structures are observed in CC processes. Differential CC dijet distributions 
are measured for the first time. Perturbative QCD calculations in NLO based on the MEPJET 
program describe the data well within errors. These calculations suggest that both quark and 
gluon contributions are needed to give a consistent description of dijet production in CC pro- 
cesses. The NC dijet distributions, measured in the same kinematic range, are well described 
by perturbative QCD predictions in NLO based on DISENT. 

The measured data sample contains events up to Q 2 ps 35 000 GeV 2 and mi 2 « 100 GeV 
and therefore probes QCD down to shortest distances. Using parton densities derived from NLO 
QCD fits to inclusive DIS cross sections, perturbative calculations based on the electroweak and 
the strong (0(a 2 )) matrix elements are found to give a consistent description of both the NC 
and CC jet cross sections at highest dijet masses and highest Q 2 . 

Comparison of the CC with the NC jet distributions confirms that at short distances the QCD 
dynamics of the hadronic final state are essentially independent of the underlying electroweak 
scattering process as is expected within the standard model. 
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0.02 
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Table 2: Rates of CC and NC events with one, two and three jets as a function of Q 2 . The 
events satisfy fi pt > 25 GeV and 0.03 < y < 0.85. The jets are reconstructed using the 
modified Durham algorithm with a fixed jet resolution parameter y cut = 0.002. The jets satisfy 
the cut 10° < 9 jet < 140°. The relative statistical errors 5 stat and relative systematic errors 5 sys 
are given in per cent. 
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Table 3: Normalized dijet cross sections as a function of y 2 , m\ 2 , z p , x p , E T j w d and 9f w d in 
CC and NC events with p l ^ pt > 25 GeV and 0.03 < y < 0.85 determined with the modified 
Durham algorithm. The events satisfy the cuts y 2 > 0.002 and 10° < 9j et < 140°. The relative 
statistical errors S stat and systematic errors 5 sys are given in per cent. 
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Table 4: Normalized dijet cross sections as a function of in CC and NC events with Q 2 > 
5000 GeV 2 , V % vt > 25 GeV and 0.03 < y < 0.85 determined with the modified Durham 
algorithm. The events satisfy the cuts y 2 > 0.002 and 10° < 9j et < 140°. The relative statistical 
errors 5 sta t and systematic errors 5 sys are given in per cent. 
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Figure 1: Feynman graphs for DIS in lowest order (a), and selected leading-order diagrams 
contributing to dijet production: Boson-Gluon-Fusion (b) and QCD-Compton scattering (c). 
The variables x and £ denote the fraction of the proton's momentum P carried by the scattered 
parton. 
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Figure 2: Uncorrected data distributions of P T , Q 2 , x and y for the selected CC and NC events. 
The observables are calculated using the hadronic final state for CC events, and the scattered 
positron for NC events. The errors are statistical only. Also shown are the corresponding 
predictions of the MC model ARIADNE 4.10 including radiative QED corrections and the HI 
detector simulation for CC processes (full line) and NC processes (dashed). 
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Figure 3: A display of two CC events. The left part shows a side view of the HI central 
and forward tracking systems surrounded by the electromagnetic and hadronic sections of the 
liquid argon calorimeter and of the lead/scintillating-fibre calorimeter. The full lines and filled 
rectangles correspond to tracks reconstructed in the tracking systems and energy depositions in 
the calorimeter, respectively. The proton beam enters from the right. The right part shows a 
view along the beam of the same events. For the upper event y 2 ~ 0.00008 and mi 2 ~ 12 GeV. 
For the the lower event y 2 ~ 0.013 and m 12 ~ 73 GeV. 
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Figure 4: Rates of CC and NC events with one, two and three jets as a function of Q 2 . The 
events satisfy p l ^ pt > 25 GeV and 0.03 < y < 0.85. The jets are reconstructed using the 
modified Durham algorithm with a fixed jet resolution parameter y cut = 0.002. The jets satisfy 
the cut 10° < 6 jet < 140°. Also shown are the predictions of the MC model ARIADNE 4.10. 
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Figure 5: Distributions of rni2, z p , x p , E T) f wd and 9j wd in CC events with p^ p > 25 
GeV and 0.03 < y < 0.85 determined with the modified Durham algorithm. The events 
satisfy the cuts y 2 > 0.002 and 10° < 6j et < 140°. Also shown are the predictions of the 
MC models ARIADNE 4.10 (full line), LEPTO 6.5.2(3 (dashed), HERWIG 5.9 (dotted) and 
RAPGAP 2.08/06 (dashed-dotted). 
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Figure 6: Distributions of y 2 , rni 2 , z p , x p , E Tj f wd and Qf w & in NC events with p^ p > 25 
GeV and 0.03 < y < 0.85 determined with the modified Durham algorithm. The events 
satisfy the cuts y 2 > 0.002 and 10° < 6j et < 140°. Also shown are the predictions of the 
MC models ARIADNE 4.10 (full line), LEPTO 6.5.2(3 (dashed), HERWIG 5.9 (dotted) and 
RAPGAP 2.08/06 (dashed-dotted). 
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Figure 7: Distributions of y 2 , mi 2 , z p , x p , E T , j w d and 6f wd in CC events with p^P > 25 GeV 
and 0.03 < y < 0.85 determined with the modified Durham algorithm. The events satisfy the 
cuts yi > 0.002 and 10° < 6j et < 140°. Also shown are perturbative QCD calculations in 
NLO obtained with MEPJET combined with a correction for hadronization effects. The shaded 
area shows the hadronization uncertainties and the renormalization scale uncertainties of the 
NLO calculations added in quadrature. In addition, the jet distributions obtained for quark- and 
gluon-induced processes are shown separately. 
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Figure 8: Distributions of y 2 , rriu, z p , x p , E T , f w d and 6j w d in NC events with p£ p > 25 GeV 
and 0.03 < y < 0.85 determined with the modified Durham algorithm. The events satisfy the 
cuts ?/2 > 0.002 and 10° < 6j et < 140°. Also shown are perturbative QCD calculations in 
NLO obtained with DISENT combined with a correction for hadronization effects. The shaded 
area shows the hadronization uncertainties and the renormalization scale uncertainties of the 
NLO calculations added in quadrature. In addition, the jet distributions obtained for quark- and 
gluon-induced processes are shown separately. 
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Figure 9: Distributions of m\i in CC events (full circles) and NC events (empty circles) 
for p l ^ pt > 25 GeV and 0.03 < y < 0.85 with and without the additional requirement 
Q 2 > 5000 GeV 2 determined with the modified Durham algorithm. The events satisfy the 
cuts y 2 > 0.002 and 10° < 9j et < 140°. Also shown are perturbative QCD calculations 
in NLO obtained with MEPJET (for CC) and DISENT (for NC) combined with a correction 
for hadronization effects. The lower histograms and the corresponding data points have been 
scaled by a factor of 1/100. The shaded area shows the hadronization uncertainties and the 
renormalization uncertainties of the NLO calculations added in quadrature. For clarity, the NC 
uncertainties are shown only. The CC uncertainties are of similar size. 
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Figure 10: The distributions of m 12 , y 2 and 6f wd in CC events (full circles) and the correspond- 
ing distribution in NC events with reweighting (empty circles) as described in the text. The 
solid histogram corresponds to the NC distributions without reweighting, which are also shown 
as data points in Figures |6| and [|. The same jet selection criteria as above are applied. 
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